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ABSTRACT 

We use observed UV through near IR spectra to examine whether SN 201 Ife can be 
understood in the framework of Branch-normal SNe la and to examine its individual 
peculiarities. As a benchmark, we use a delayed-detonation model with a progenitor 
metallicity of Zq/20. We study the sensitivity of features to variations in progenitor 
metallicity, the outer density profile, and the distribution of radioactive nickel. The 
effect of metallicity variations in the progenitor have a relatively small effect on the 
synthetic spectra. We also find that the abundance stratification of SN 201 Ife resem¬ 
bles closely that of a delayed detonation model with a transition density that has been 
fit to other Branch-normal Type la supernovae. At early times, the model photosphere 
is formed in material with velocities that are too high, indicating that the photosphere 
recedes too slowly or that SN 20life has a lower specific energy in the outer « 0.1 Mq 
than does the model. We discuss several explanations for the discrepancies. Finally, 
we examine variations in both the spectral energy distribution and in the colors due 
to variations in the progenitor metallicity, which suggests that colors are only weak in¬ 
dicators for the progenitor metallicity, in the particular explosion model that we have 
studied. We do find that the flux in the U band is significantly higher at maximum 
light in the solar metallicity model than in the lower metallicity model and the lower 
metallicity model much better matches the observed spectrum. 

Key words: radiative transfer - supernovae: general - supernova: individual: 
SN 2011fe. 


1 INTRODUCTION 

2012 

jVinko et al.||2012| [Munari et al.||2013| 

Pereira et al. 


2013 

), and in the IR (Matheson et al. [2012). 

Pereira et al. 


Supernova PTFllkly/2011fe (henceforth SN 2011fe) was 
discovered by the Palomar Transient Factory on Aug 23, 
2011 in MlOl only hours after explosion ( Nugent et al.|20lT ). 
This nearby, early discovered object has been extremely well 
observed in all bands. The Carnegie Supernova Project ob¬ 


tained excellent spectroscopy in the infrared (Hsiao et al. 
2013). Photometry has been obtained by SWIFT in the 
UV ( Brown et al.|2012 ), in the optical ( [Richmond Smit h| 


(2013) presented a detailed comparison of the photometric 


observations including a well calibrated spectrophotomet- 
ric time series. SN 201 Ife has been observed in the radio 


(Chomiuk et al. 2012), in gamma-rays (Isern et al. 2013) 


and in the X-ray (Liu et al. 2012). From the spectra and 


photometry SN 201 Ife is about as ordinary a SN la as there 
could be. Due to the early discovery and close proximity, 
several groups have drawn conclusions about the environ- 
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merit of SN 201 Ife. Using non-detections and the very early 
observed points [Bloom et al. (2012) were able to constrain 
the radius of the primary star, concluding that it must be a 
compact object (white dwarf or neutron star). Using archival 
HST images [IT et al.[ ( [2011[ ) were able to rule out luminous 
red giant and helium star companions. Using radio data, 
constraints have been placed on the progenitor environment 
which have been interpreted as ruling out the single degen¬ 


erate scenario (Chomiuk et al. 2012); however, the winds 


blown during the progenitor formation could naturally pro¬ 
duce a low density environment that do not necessarily re¬ 
quire a degenerate companion (see also Horesh et al.|20l'2 ). 

Some spectral modeling of the optical spectra of 
SN 201 Ife has been presented ( [Ropke et al. 2012 Dessart 


et al. 2014). Here, we want to compare and contrast 


SN 201 Ife to a model which reproduces the spectra of a 
Branch-normal supernova ( Sandage et al.[[T996 ). We use a 
model because this allows us to translate spectral similari¬ 
ties and differences into physical space, for example, to de¬ 
termine the relevant mass layer involved in the differences. 
We have specifically chosen a generic DD model for Branch- 
normal supernova and not tuned the model for this specific 
supernova. The goal of this study is to show that SN 20life 
is close to a Branch-normal supernova and to evaluate dif¬ 
ferences and discuss possible physical causes. 

This study is based on the combined UV and opti¬ 
cal HST and IR Gemini spectra. The infrared spectra that 
we study have been already were presented in | Hsiao et al.| 
(2013). The UV+optical observations were also presented 
and studied in previous work ( Foley et al.|2012a]b Foley & 


Kirshner]|2013 Mazzali et al.[|2014 ). Here, we present a de¬ 

tailed comparison quantitative synthetic spectroscopy with 
a number of epochs of SN 201 Ife, all of which include cov¬ 
erage from the UV through the IR. 


2 SPECTRAL CALCULATIONS AND 
EXPLOSION MODEL 

We performed spectral calculations using the multi¬ 
purpose stellar atmospheres program PHOENIX/ID version 
16 (Hauschildt V Baron 1999[ [Baron V Hauschildt 1998 


Hauschildt, Baron fc Allard] 1 1 997[ [Hauschildt et^ff 1997 

1996). Version 16 incorporates many changes over previ¬ 


ous versions used for supernova modeling (Baron, Branch 
fc Hauschildt[2007 Baron et al.|2006" ) including many more 
species in the equation of state (83 versus 40), twice as many 
atomic lines, many more species treated in full NLTE, and 
an improved equation of state. PHOENIX/ID solves the radia¬ 
tive transfer equation along characteristic rays in spherical 
symmetry including all special relativistic effects. The non- 
LTE (NLTE) rate equations for many ionization states are 
solved, including the effects of ionization due to non-thermal 
electrons from the 7-rays produced by the radiative decay of 
^®Ni, which is synthesized in the supernova explosion. The 
atoms and ions calculated in NLTE are He I-H, C I-HI, 
O I-HI, Ne I, Na I-H, Mg I-HI, Si I-HI, S I-HI, Ca H, 
Ti H, Cr I-HI, Mn I-HI, Ee I-HI, Co I-HI, and Ni I-HI. 
These are all the elements whose features make important 
contributions to the observed spectral features in SNe la. 

Each model atom includes primary NLTE transitions, 
which are used to calculate the level populations and opacity. 


and weaker secondary LTE transitions which are included in 
the opacity and implicitly affect the rate equations via their 


effect on the solution to the transport equation (Hauschildt 


& Baron 1999). In addition to the NLTE transitions, all 


other LTE line opacities for atomic species not treated in 
NLTE are treated with the equivalent two-level atom source 


function, using a thermalization parameter, a = O.IO (Baron 
et al.|I99'6 ). The atmospheres are iterated to energy balance 
in the co-moving frame; while we neglect the explicit effects 
of time dependence in the radiation transport equation, we 
do implicitly include these effects, via explicitly including 
p dV work and the rate of gamma-ray deposition in the gen¬ 
eralized equation of radiative equilibrium and in the rate 
equations for the NLTE populations. 

The outer boundary condition is the total bolometric lu¬ 
minosity in the observer’s frame. The inner boundary condi¬ 
tion is that the fiux at the innermost zone {v — 700 km s“^) 
is given by the diffusion equation. Converged models re¬ 
quired 256 optical depth points to correctly obtain the Si H 
A6355 profile. 

Eor our analysis, we use a delayed-detonation model 
which reproduces the light curves and spectra for Branch- 


normal supernovae (Hoefiich, Khokhlov & Wheeler| 

1995 

Hoefiich et al. 

|2002 

Hoefiich 2006 Gerardy et al. [2003). 

Gerardy et al. 

to 

0 

0 

CO 

) compared a very similar model to 


The models start from a C/0 white dwarf taken from the 
core of an evolved SM© main sequence star. Through ac¬ 
cretion, this core approaches the Chandrasekhar limit. An 
explosion begins spontaneously when the core has a c entral 
density of 2.0x10^ gcm“^ and a mass close to 1.37 M© (Hoe- 
fiich||20Q2 ). The transition from defiagration to detonation 
is triggered at a density of 2.3 x 10^ g cm“^. We considered 
two modes for the delayed-detonation transition: one with a 
direct transition during the defiagration phase, the other af¬ 
ter a mild pulsation which formed an envelope of 10“^ M©. 
The resulting density structures are shown in Eig. We 
considered initial metallicities Z of Zq and Z©/20. Here, Z 
is defined as the iron abundance relative to solar. We take 
into account the smaller variation of the elements such as 


Ne and O compared to the Ee-group elements (Argast et al 


2001). The former dominates the metallicity effect on nu¬ 


clear burning, whereas the latter sets the fioor for Ee-group 


elements (Hoefiich, Wheeler & Thielemann 
et al.|2000| |. 


1998 


Hoefiich 


3 RESULTS 

3.1 General Properties 

We calculated both a classical delayed detonation (DD) 
model and a pulsating delayed detonation (FDD) model mo¬ 
tivated by the report of Parrent et al. ( 20I2| of two distinct 
high and low velocity components in the Si, O, and C lines. 
In Eigure[^ we show the comparison between the DD model, 
the FDD model, and the observations at day 2. While nei¬ 
ther model is strongly preferred over the other, the features 
in the PDD model are more washed out, a trend which con¬ 
tinues into later epochs and it indicates that the dense shell 
is certainly not desirable, thus we choose to use the standard 
DD model for this study. Our PDD model and DDT model 
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Models for SN 201 Ife 3 




Figure 1. The density profile of the prompt DDT and the FDD 
that were compared. The value of po at maximum light is about 
2 X 10“^^ g cm“^. 



A (A) 


Figure 2. FDD vs DD at Day 02. The spectra of the FDD model 
(dashed line) and the DD model (dot-dashed line) are compared 
to the observations (solid line). 


differ somewhat from those of Dessart et al. (2014), but we 
find that the variation in the colors, specifically the flux in 
the D—band is due to metallicity of the progenitor. 

We find a preference between models with solar metal¬ 
licity and Zq/20. There are not strong differences at the ear¬ 
liest times, but rather at maximum light. The solar metal¬ 
licity progenitor produces somewhat redder spectra in the 
bluest bands and we adopted a progenitor metallicity with 
Z0/2O as the hducial model. We discuss the differences in 

§in] 


3.2 Benchmarking SN 201 Ife 

Having probed the primordial metallicity and explosion 
class, we use a DD model whose calculated light curve gives 
results in line with the observed class of Branch-normal su¬ 
pernovae as a benchmark for SN 201 Ife. With time, spec¬ 
tra reveal increasingly deeper layers. The spectral sequence 
provides a key probe of the layers for similarities and differ¬ 
ences between the model and the observations. Therefore, 
throughout our discussion we will identify the layers by their 
mass coordinate measured from the outside — more pre¬ 
cisely, we report the location of the line forming region of 
the Si II A6355 line. The velocity of both the observed and 
synthetic line is measured the same way, by determining the 
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Figure 3. The bolometric light curve from these calculations 
compared to the bolometric light curve inferred by|Fereira et al.| 
( | 201 ^ . 


wavelength of the blueshifted minimum. Note that the UV 
spectra always probe the very outer layers (see Fig. 11 in 
Hoehich, Khokhlov Wheeler [1 995 ). 

We consider epochs up to maximum light because the 
corresponding “photospheric regions” undergo partial burn¬ 
ing. Details of the spectral features are rather sensitive to 
small variations in brightness and non-thermal excitation 
( Baron et al.||2012 ). 

3.2.0.1 Aug 25/Day 2 The earliest spectrum of 
SN 201 Ife was obtained on Aug. 25th, about 2 days af¬ 
ter the explosion, which probes the outer 5 x 10~^Mq. The 
spectral features are quite sensitive to the temperature and 
excitation, we show the NLTE spectra in Figure We cal¬ 
culated an LTE model with the same parameters and the 
LTE model is much too bright in the IR and the features 
are weak, compared to the NLTE model. The hgures show 
that almost ah the lines present in the observed spectrum 
are also in the synthetic spectrum. While many of the line 
strengths are well reproduced by the model, the Si II A6355 
line is much too weak and the Si II A5970 line is about the 
right strength, but is too fast. The LTE model over-predicts 
the IR hux and the IR features are weak, whereas the NLTE 
model roughly gets the IR continuum correct, but the IR fea¬ 
tures are too strong and broad. A common feature of both 
the LTE and NLTE spectra is a too large Doppler shift of 
ah the absorption features of intermediate mass elements 
both in the optical and IR, including the Ca II IR-triplet 
and the O I A7773.4 line. As mentioned above, details of the 
line-features are very sensitive to excitation and tempera¬ 
ture effects. In particular, a higher Doppler-shift may be 
produced by an excessive emission component seen for ex¬ 
ample in Mg II AI0926 or O I A7773.4. However, even in case 
of LTE, the trend is conhrmed. The smaller Doppler shifts 
in SN 201 Ife compared to the delayed detonation model for 
a Branch-normal supernova may be attributed to a lower 
specihc energy or photosphere that recedes more quickly. 

3.2.0.2 Aug 28/Day 5 Eigure|^ shows the comparison 
to the combined i75'T+optical+IR spectrum obtained on 
MJD 55801.12 which probes the outer 3 x 10~^Mq. The 
overall spectral shape from UV through optical is quite well 
reproduced. The far UV is a bit high, but the quality of 
the ht from 0.2 — 2.5 /im is quite good, the Si H A6355 
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Figure 4. NLTE synthetic spectrum on day 02 (dashed line) is 
compared to the optical spectrum obtained at Lick on Aug 25, 
2011 and the IR spectrum obtained at Gemini North on Aug 25, 
2011 (solid line). 



A (A) 


Figure 6. NLTE synthetic spectrum on day 08 (dashed line) is 
compared to the UV/optical spectrum obtained by HST on on 
Aug 31, 2011 (MJD 55804.25) and the IR spectrum obtained at 
Gemini North on Aug 31, 2011 (solid line). 



5000 10000 15000 

A (A) 

Figure 5. NLTE synthetic spectrum on day 05 (dashed line) is 
compared to the UV/optical spectrum obtained by HST on on 
Aug 28, 2011 (MJD 55801.12) and the IR spectrum obtained at 
Gemini North on Aug 28, 2011 (solid line). 


line is about the right strength, the Si II A5970 line is too 
strong and fast by about 4,000 km s“^, as is the Ca IR 
triplet. Actually, the the Si II A5970 line is low by about 20 
percent, but it does not recover to the blue as the observed 
spectrum does, due to blending with other lines. In general 
the spectrum is somewhat too fast. 


Al 1714.87 feature is prominent in the synthetic spectrum, 
but much weaker in the observed spectrum. 

3.2.0.4 Sep 10/Day 18 In Fig. we show the spec¬ 
trum obtained on MJD 55814.43 which which is close to 
maximum light (Nugent et al. 2011| [Pereira et al.||2013 ) 
which probes the outer O.bM©. For reference, [Pereira et al.j 
(2013) find on MJD 55814.51. The observed spectrum 


of SN 201 Ife closely resembles the synthetic spectrum for 
a Branch-normal supernova both with respect to the line 
strength and the Doppler shifts of individual lines. Overall 
the Doppler shifts are consistent, but the model is slightly 
too cool, resulting in a reduced flux in the blue. This could 
be due either to strong line blanketing in the blue, which 
pushes flux to the red or because the opacities due to the 
iron group elements decrease in the blue as the temperature 
decreases. The complex formation of the observed spectrum 


was studied in Bongard et al. (2008). The Si II A6355 line 


is now in good agreement in both absorption and emission 
strength, indicating that the model at this phase has about 
the correct velocity extension. The Si IIA5970 line now is also 
well fit in both absorption and emission (modulo the too low 
pseudo continuum). The Ca II IR triplet is too weak. The 
Mg II A9226 and Mg II A10926 features are too strong, as 
are most of the rest of the features further to the red in the 
IR. 


3.2.0.3 Aug 31/Day 8 Figurej^shows the comparison 
to the combined i75'T+optical+IR spectrum obtained on 
MJD 55804.25 which probes the outer 5 x 10“^M©. Not 
only is the overall spectral shape from UV through optical 
quite well reproduced, but also the relative strength of the 
far UV closely matches the observations. However, the iron 
and silicon feature around 4000 A is too weak. The strength 
of Ca II H+K is well reproduced. The calcium IR triplet is 
now too narrow, and not strong enough in absorption. The 
Si II A6355 line is well fit in absorption strength, but it is 
too fast by about 4,000 km s“^ and the emission peak is a 
about 20% too high, indicating that the model at this phase 
is too extended. The Si II A5970 line is now also well fit in 
absorption. The Mg II A9226 line is much too strong. The 
Mg II A10926 feature is too strong as are most of the rest of 
the features further to the red in the IR. The O I A7773.4 
line is prominent, but too strong and too broad. The Si II 


4 DISCUSSIONS AND CONCLUSION 

Figure shows the optical spectra of the Aug 25 spectrum 
we have modeled, along with the synthetic spectra shifted 
by a “redshift” z = 0.015. Similar experiments give “redshift 
values of 0.015,0.015,0.010, for Aug 28, Aug 31, and Sep 
10, respectively. This shows that the recession of the pseudo 
photosphere in the hydro model with time is slower than in 
the observations. We can estimate the amount of the shift 
by comparing the velocity of the absorption minimum of 
Si II A6355 line at each epoch. Table shows the measured 
minimum velocity for each epoch. Since the velocity is close 
at maximum light we will consider the “outer part” of the 
model to be at velocities v ^ 11,300 km s“^. 

We can reduce the effective velocity extension by steep- 
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Models for SN 201 Ife 5 



Figure 7. NLTE synthetic spectrum on day 18 (dashed line) is 
compared to the UV/optical spectrum obtained by HST on on 
Sep 10, 2011 (MJD 55814.43) and the IR spectrum obtained at 
Gemini North on Sep 10, 2011 (solid line). 



Figure 9. The original density profile of the hydro model is com¬ 
pared to a modified profile where the density is forced to fol¬ 
low a powerlaw p oc (r’/r’cut)~^, for velocities v > Vcut, where 
Vcut = 12, 000 km s“^, and n = 12. 


Si II A6355 Velocities 


Epoch (Day) 
(since explosion) 

Fsiii 

(observed) 

Fsiii 

(synthetic) 

2 

16,300 

20,700 

5 

13,000 

17,300 

8 

11,300 

16,000 

18 

9,700 

10,100 


Table 1. Velocities of the absorption minimum of Si II A6355 
feature at each epoch in km s“^. 


ening the density profile beyond the photosphere. This can 
be produced by an outwardly increasing specific energy with 
mass element. While steeper density profiles violate energy 
conservation from nuclear burning and make the models hy- 
drodynamically inconsistent, we perform an empirical ex¬ 
ercise to examine its effect on the synthetic spectra, being 
careful not to generalize the results too much. Figurej^shows 
the modified density profile, obtained by forcing the density 
to follow a power-law p oc {v/vcut)~^, for velocities v > Vcut , 
where we took I’cut = 12,000 km s“^, and n = 12 (Branch 


Branch, Dang fc Baron||2009| 


et al. 


Doull Baron[|2011 ). Figure [lO] shows a comparison of the 



A (A) 


Figure 8. Optical blow up of Aug 25 spectrum (solid line) shown 
in Fig where the synthetic spectra are shown with (dot-dashed 
line) and without (dashed line) a “redshift” of 2 ; = 0.015. 



A (A) 


Figure 10. NLTE spectra day 05 original hydro model (dashed 
line) compared to model hydro model with density profile p oc 
v~^‘^ for V > 12,000 km s“^ (dot-dashed line). The observed 
spectrum is shown with a solid line. 


Aug 28/Day 5 spectrum with and without the density pro¬ 
file modification. The effects of the density modification are 
largest in the red and smaller in the blue. The Si IIA6355 is 
somewhat slower and the line profile is clearly narrower, but 
the agreement is not significantly better. However, in the 
red, the shape of the O I A7773.4 line is much better and 
that continues on to the Ca II IR triplet, the Mg II A9226 
feature, and the Mg II A10926 feature. In the modified den¬ 
sity structure the Si II Al 1714.87 becomes more pronounced. 
Thus, while steepening the density profile does not signifi¬ 
cantly improve the agreement, it does have some benefit in 
moving the absorptions of the redder lines to the correct 
velocity. 

A full study of the possible reasons for lower observed 
photospheric velocities is beyond the scope of the current 
study. Some possible reasons for lower velocities can be un¬ 
derstood within the framework of spherical DD models. For 
a wide range of model parameters the overall density struc¬ 
tures are very similar (Hoeflich, Muller & Khokhlov 1993). A 
possible reason for the lower observed photospheric velocity 
may be a faster receding photosphere in mass. To first order 
for one to two weeks past maximum light in Branch-normal 
SNe la, the opacities remain high and the photosphere re¬ 
cedes in mass due to geometrical dilution of the expand- 
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ing envelope ( Hoeflich, Muller fc Khokhlov||1993 ). In order 
to increase the recession rate of the pseudo-photosphere we 
need a lower kinetic energy. In classical DD models, most of 
the matter undergoes burning. Possible ways to reduce the 
kinetic energy include: higher binding energy of the white 
dwarf and thus, higher central densities; less nuclear energy 
production due to a smaller C/0 ratio, by, for example, a 
larger main sequence mass; less heating in the outer layers 
by ^®Ni by shifting its distribution to the central layers. We 
are not signaling out one particular cause and to do so is be¬ 
yond the scope of this work. We note that Baron et al. (2012) 
suggested that a higher value of C/0 could explain the high 
velocities seen in SN 2001ay, whereas Ohlmann et al. (2014) 
found that varying the C/0 ratio alone in a fixed white 
dwarf did not have dramatic effects. 

PDD models with a large amount of unburned material 
will reduce the explosion energy proportional to the amount 
of unburned material, and may lower the velocities of el¬ 
ements of incomplete burning. However, strong pulsations 
produce an outer, high velocity layer as observed in, for ex¬ 
ample, SN 1990N and similar SNe la ( Quimby et al.||20Q6 ). 
An alternative explanation is significantly lower mass pro¬ 
genitor than the Chandrasekhar mass which, would be less 
luminous and inconsistent with the light curve. 


4.1 Progenitor Metallicity 

Figure pT] compares the synthetic spectra of the fiducial pro¬ 
genitor model with Z — Z 0 / 2 O to that of progenitor with 
solar metallicity. The other explosion parameters were the 
same for both models. By eye, prior to maximum light, there 
is not a clear preference for either model, which is somewhat 
surprising since naively one expects progenitor metallicity to 
play a role at particularly early times. At maximum light a 
strong discrepancy appears, where the solar metallicity pro¬ 
genitor appears too blue or more accurately, the flux is much 
higher in the f/—band. To attempt to quantify this we calcu¬ 
lated synthetic photometry on both the observations and the 
models and compared the colors. This procedure is subject 
to systematic errors since it assumes that the relative flux of 
the observations is very accurate, so that the synthetic colors 
calculated from the observed spectra are meaningful. Given 
this caveat, we did not find a clear preference for either pro¬ 
genitor metallicity, based solely on colors. We examined the 
25 filters from HST and the Johnson set UBVIRJHKs and 
looked at all combinations of neighboring blue filter - red fil¬ 
ter. The lack of clear trends can be seen in the bluest filters 
F220W - F250W, the models are too red at day 02, and too 
blue at days 08 and 18, in the F250W - F330W filters, the 
models are significantly too blue at day 02, but then only 
slightly too blue at days 08 and 18. The discrepancy is far 
smaller in the F330W- F344N and F344N - F435W colors, 
with the same trend and then changes in the F435W - F475 
color where the models are too blue at day 02, too red at day 
08 and then the solar model is slightly too blue and the sub¬ 
solar model slightly too red at day 18. Figure shows the 
differences in the twenty-three different colors at the three 
different epochs studied. The trend for the solar metallicity 
model to become much redder in the bluest bands at maxi¬ 
mum light is evident, but it is difficult to discern an overall 
pattern in the colors. In the bluest filters the solar metallic¬ 
ity model is much redder at maximum light, until we reach 



A (A) 


Figure 11. The synthetic spectra of the Z = ZqI2{) (dashed line) 
and the Z = ZqI2{) (dot-dashed line) progenitor model compared 
to observations (solid line) for days Aug 25 (top panel), Aug 31 
(middle panel), and Sep 10 (bottom panel). 


the F435W -F475W color, where the solar metallicity model 
suddenly becomes bluer than the lower metallicity model, 
red ward there is no strong trend in the colors. Prior to max¬ 
imum light, there is really no strong trend in the colors at 
any wavelength. The largest discrepancy is in fact evident 
from Figure m where the solar metallicity model is clearly 
much brighter in the f/—band. In fact, the solar metallicity 
model is nearly 0.5 mag brighter than the low metallicity 
model in f/, even though both models are blue mU — B and 
the solar metallicity model is only about 0.1 mag bluer in 
U — B than the low metallicity model. Thus, one really needs 
spectra in the restframe f/—band to try to get a handle on 
progenitor metallicity. 

Overall, the results make sense, in the bluest colors, 
particularly at maximum light, the solar metallicity model is 
much redder than the lower metallicity model due to higher 
opacities in the UV, but clear diagnostics are not so evident. 
Our results are in good agreement with those of Hoeflich 


and collaborators (Hoeflich, Wheeler & Thielemann 1998 
Hoeflich et al.|2000 ). Mazzali et al.| ( |20T4 ) found evidence for 

a progenitor metallicity of ^o/2, which is not inconsistent 
with our results. 

Figure shows that the C H A6580 appears with 
roughly the strength and shape as does the observed fea¬ 
ture in the solar metallicity model at on Aug 25, while not 
at the correct Doppler shift. The velocities of the feature are 
(15, 000,18, 000,19, 000) km s“^ for the observations, the so¬ 
lar metallicity progenitor, and the low metallicity progeni¬ 
tor, respectively. We have confirmed the identification of the 
line as due to C H in the models by rerunning the spectra 
with the C H line opacity set to zero. This gives credence 


to claims of C H in other SNe la (see [Thomas et al.||2011 

Parrent et al. 

|2011| [Thomas et al.|2007[[Folatelli et al.|2012 

Parrent||2015| 

, and references therein). 


in terms of the variation in the progenitor evolution. At the 
end of the central He-burning and therefore, in mixtures 
which are helium poor, ^^C(a,7)^®0 results in lower C/0 
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Figure 12. The difference in 23 synthetic colors for z = Zq 
the Z = progenitor model compared for days Aug 

25, Aug 31, and Sep 10. The colors are plotted from bluest 
to reddest filters and the difference plotted is, for example, 
(B — V)zq — (B — V)zQf20- The filters used are in order 
from blue to red, HST/ACS-HRC: F220W, F250W, F330W, 
F344N, F435W, F475W, F502N, F550M, F555W, F606W, 
F625W, F658N, F660N, F775W, F814W, F850LP, F892N; UB- 
VR1+2MASS: U,B,V,R,I,J,H,K. We did not calculate colors 
across the filter sets, so the F892N-U color is not plotted. 



A (A) 


Figure 13. The synthetic spectra of the Z = Zq (dashed line) 
and the Z = Zq/20 (dot-dashed line) progenitor model compared 
to the observation (solid line) on Aug 25. While both models show 
a line due to C II A6580, the solar metallicity model line is about 
the right shape and strength as the observed feature. However the 
velocity minima are at 15,000 km s“^ in the observations, 18,000 
km s“^ in the solar metallicity model and 19,000 km s“^ in the 
model with Z = ZqI2{). 


ratios in the convective core. With lower metallicity, the size 
of the convective core increases due to lower opacities. Shell 
burning will produce material with C/O ~ 1. As a result 
of less carbon, the explosion energy for our progenitors de¬ 
creases with Z. With higher explosion energy, we expand 
faster and for a given DDT, the very outer layers are burned 
at higher density reducing the amount of material in the un¬ 
burned layers and resulting in a higher expansion velocities 


(Hoeflich, Wheeler & Thielemann|1998 Hoeflich et al.|2000 

2oTor 



Figure 14. The deposition function at day 02, using the model 
nickel distribution and the one obtained by using the prescription 
of [Ph^ ( [m^ . 


4.2 Enhanced Nickel Mixing? 


We examine a suggestion of Piro (2012) who by studying the 
early rise of the light curve concluded that ^®Ni was required 
in the outer 0.1 < M< 1 x 10“^ of the supernova. He fur¬ 
ther concluded that this requirement places constraints on 
the explosion model, somewhat favoring models where the 
detonation begins in the outer parts of the star. Piro’s [20T^ 
study is based on the diffusion time scales for the rise time. 
In fact, geometrical dilution will be responsible for the re¬ 
ceding of the photosphere and adiabatic expansion of the 
corresponding layers must be taken into account. Figure 
shows the 7-ray deposition function (the 7-ray deposition 
function is the ratio of the rate of absorption of 7-rays J 
to the instantaneous decay rate of ^®Ni, see, for example, 
Sutherland fc Wheeler || 1984 Jeffery || 1998 ) for our standard 
case, compared to that obtained using the prescription of 
Piro (2012). We chose to make the mass fraction of radioac¬ 


tive nickel constant at = 0-03 above a velocity of 

15,500 km/s, which corresponds to a mass of 0.1 M© masses 
from the outer boundary. This is about in the middle of the 


parameter choices outlined in Piro (2012). In our model the 
bottom layer where the ^^Ni is added occurs at a velocity of 
V ~ 15,500 km s“^. Thus, the photosphere in our model is 
signihcantly above this modihed heating region. Neverthe¬ 
less, we hnd negligible changes due to the added heating. 
The temperature of our models is a little higher in the outer 
parts, but not enough higher in the low density environment 
to signihcantly alter the observed spectra. Figure [^explains 
this result: the 7-ray deposition function 0dep is increased by 
the extra nickel only up to the level of ~ 1 x 10~^, but only 
at the highest velocities and not near the photosphere. The 
increase in deposition is not enough to dramatically change 
the outer energy deposition rate at Day 02. 


4.3 Comparison to other work 


Light curves and spectra of SN 201 Ife at these epochs have 


been studied by 

Ropke et al. ( 

2012), Foley (2012), Foley & 

Kirshner 

(2013), 

Dessart et al. 

(2014), Mazzali et al. ( 

2014), 

and Graham et 

al. (2015). Ropke et al. (2012) compared 


the hdelity of a delayed detonation Chandrasekhar mass 
model and violent merger model to the observations of the 
Nearby Supernova Factory and found that both models had 
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their shortcomings, but found no strong preference for ei¬ 
ther model. Mazzali et ah (2014) performed abundance to¬ 
mography on a W7 model and a delayed detonation model 
and found that the delayed detonation model was some¬ 
what preferable, in addition to finding a primordial metal- 
licity that was about half solar. Dessart et al. (2014) studied 
PDD models in general, and found a preference for a PDD 
model to match the B-band light curve of SNe 201 Ife and the 
spectral evolution. However, their PDD models have much 
weaker shells then the model we considered (compare their 
Figure 10 to Figure [^. They also find very washed out fea¬ 
tures at early epochs (see their Figure 16). [Foley Kirshner 


(|2013| ) used empirical comparisons with the models of Lentz 


et al.| ( [20QQ ) to find a low metallicity for SN 201 Ife of about 


Zq/4 and a ratio between the metalicity of SN 2011by and 
SN 2011fe of about 30. Graham et al. (2015) found the that 
UV flux in the assumed supersolar SN 201 Iby is significantly 
lower than that in SN 201 Ife, in rough agreement with the 
W7-based models of Lentz et al. (2000), but in contradis¬ 
tinction to the DD models used in this work. Thus, overall 
our results are generally compatible with those of previous 
work. 


5 CONCLUSIONS 


We have compared HST and ground based spectra of 
SN 20life up to maximum light to detailed NLTE synthetic 
spectra of a Branch-normal hydrodynamic model. Progeni¬ 
tor models with metallicities of of Ze and Z0/2O qualita¬ 
tively show overall reasonable results after day 2 based on 
optical spectra. 

We explored the possibility that the high velocity fea¬ 


ture reported by Parrent et al. (2012) may be produced by 


a shell as a result of a “low amplitude pulsation”. Though 
a shell may produce a high velocity feature, it also would 
produce high fluxes in the UV and U wavelengths inconsis¬ 
tent with the observations as well as the very high velocity 
photosphere at early times, which leads to washed out fea¬ 
tures (§|^and Dessart et al.|[2M4 ). Branch normal models 
show features at the corresponding wavelengths produced 
by iron-group elements ( Gerardy et al.||2003 ). 

Nearly all the features in the observed spectra are seen 
in the synthetic spectra in approximately the correct range 
of the velocity, that is, the elemental stratification in the 
model also appears in the observations. However, as dis¬ 
cussed in §3.2[ the photosphere in the model is formed at 
too high velocities at early times corresponding to the out¬ 
ermost 0.05 to 0.1 M0. In principle, this can be corrected by 
a smaller kinetic energy of the explosion, e.g., by higher cen¬ 
tral densities and, thus, higher binding energy, or high G/0 
ratios in the progenitor. Alternatively, the photosphere of 
SN 201 Ife may recede faster due to lower excitation of ions 
(from gamma-rays or the radiation field) and/or variations 
in the density structure. In SN 201 Ife, the second explana¬ 
tion is favored because, by maximum light, the photospheric 
velocities of the model and SN 201 Ife are in close agreement 
(see Table [^. 

There are a multiple ways that the adopted model could 
be adjusted to reproduce the observed discrepancies, includ¬ 
ing variations in the chemical structure, rotation of the pro¬ 
genitor, or pulsations prior and during the explosion (al¬ 


though this would have to done in accord with constraints 
from observed spectra). Future studies will examine these 
effects in detail. 
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